Uncovering the formation of ultra-compact dwarf galaxies by 
multivariate statistical analysis 
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ABSTRACT 

We present a statistical analysis of the properties of a large sample of dynamically hot old 
stellar systems, from globular clusters to giant ellipticals, which was performed in order to inves- 
tigate the origin of ultra-compact dwarf galaxies. The data were mostly drawn from Forbes et 
al. (2008). We recalculated some of the effective radii, computed mean surface brightnesses and 
mass-to-light-ratios, estimated ages and metallicities. We completed the sample with globular 
clusters of M31. We used a multivariate statistical technique (K-Means clustering), together with 
a new algorithm (Gap Statistics) for finding the optimum number of homogeneous sub-groups 
in the sample, using a total of six parameters (absolute magnitude, effective radius, virial mass- 
to-light ratio, stellar mass-to-light ratio and metallicity) . We found six groups. FKl and FK5 
are composed of high- and low-mass elliptical galaxies respectively. FK3 and FK6 are composed 
of high-metallicity and low-metallicity objects, respectively, and both include globular clusters 
and ultra-compact dwarf galaxies. Two very small groups, FK2 and FK4, are composed of Local 
Group dwarf spheroidals. Our groups differ in their mean masses and virial mass-to-light ratios. 
The relations between these two parameters are also different for the various groups. The prob- 
ability density distributions of metallicity for the four groups of galaxies is similar to that of the 
globular clusters and UCDs. The brightest low-metallicity globular clusters and ultra-compact 
dwarf galaxies tend to follow the mass-metallicity relation like elliptical galaxies. The objects of 
FK3 are more metal-rich per unit effective luminosity density than high-mass ellipticals. 

Subject headings: galaxies; giants and dwarfs - methods: data analysis - methods: statistical 
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1. Introduction 

The variety of astrophysical structures in the 
Universe, from galaxies and galaxy clusters to 
stellar remnants, is well described by essential 
physical principles ( Padmanabhan 2000f) . How- 



ever, their origin, in particular that of globu- 
lar clusters (hereafter GCs), is still a matter of 
debate. GCs are an intermediate cell of struc- 
ture between stars and galaxies, and their forma- 



ing the Universe ( 


Peebles 19691: Ashman & Zeoi 


1992; Harris et al. 


19951: Cote et al. 1998). If we 



consider star clusters as a single class of astro- 
physical objects, there are many well-known and 
poorly understood phenomena, which make their 
origin enigmatic. For example, the color bimodal- 
ity of GC systems existing in most galaxies, the 
absence of a clear mass-metallicity relation for 
the population of red GCs, the correlation be- 
tween color and integrate d magnitude among th e 
brighter metal-poor GCs (IStrader fc SmithlbOOSi) . 



the differences in luminosity functions and surface 
density pr ofiles between young and old cluster sys- 
tems fe.g. lBrodie fc Stradeij[200l iLee et al.ll2010l 



and references therein). Spherical stellar sys- 
tems, whether globular clusters, elliptical galax- 
ies, or substructu res of spiral ga laxies, are con- 



sidered virialized (|Antonovlll973l ). The origin of 



GCs is ultimately linked to the evolution of larger 
pressure-supported st ructures within th e cosmo- 
logical hierarchy (e.g. iHwang et al. 20081 and ref- 
erences therein). 

In the last decade a new type of astronomical 



physicists (Hilker et al.ll999; 


PhilliDDS et al.|l2001 


iDrinkwater et al. 200ol 2003 


; Mieske et al.ll2006 



while making a spectroscopic survey in the Fornax 
cluster. These objects called ultra-compact dwarf 
galaxies (UCDs), dwarf globular transition ob- 
jects or sometimes intermediate massive objects, 
are different from the classical globular clusters 
or dwarf elliptical galaxies in terms of their radii, 
relaxation time and V-band mass-to-light ratios. 
They are more massive, more luminous, and have 
higher mass-to-light ratio than globular clusters, 
but are fainter and more compact than dwarf el- 
liptical galaxies. 

Several formation scenario have been pro- 
posed for understanding their physical properties. 



Kroupal (|l998[ ) and iFellhauer fc Kroupal (|2002l ) 
suggested that UCDs are the results of merger 
of many young star cluster s formed in g alaxy - 



galaxy encounters whereas iMieske et al.l (|2002r ) 
suggested that they are the luminous extension 
of massive GCs. The formation of UCDs from 
the mass threshing of the envelope s of nucleated 



galaxies has als o been suggested (IBassino et al 
1994; Zinncckcr^^talJ E988; Bc kki et al.l l2003t 
fCoerdt et al.l 120081) while along another line of 
thought UCDs are consi dered fundamenta l build - 
ing blocks of galaxies ( Drinkwater et al\ l2004[ ). 
Special efforts were made to unite old dynam- 
ically hot stellar systems, from GCs, UCDs 
and dwa rf spheroidals (dSphs) to giant elliptical 
galaxies (IZaritskv et al.ll2006t iForbes et all (20081: 
Dabringh ausen et al."2008 ) to re veal the nature of 



UCDs. (Mieske & KroupJ lioOsh have studied the 
internal dynamics of a large sample of UCDs in 
Fornax. They argue that UCDs are dynamically 
unrelaxed and dynamical evolution has probably 
not influenced their present dynamical M/L ratio. 

All these findings originated while studying 
two-point correlations between different projec- 
tions of the fundamental plane of galaxies (FP) 
defined by velocity dispersion, size (or effective ra- 
dius) and surface brightness (or mass density). For 
example, relations were found between size and lu- 
minosity, mass and metallicity, mass-to-light ratio 
and dynamical mass, luminosity and velocity dis- 
persion etc. Considering two parameters at a time 
means disregarding the combining effects of others 
which in turn are responsible for losing significant 
information. 

For a unique and robust theory of the forma- 
tion of UCDs, a multivariate approach is more ap- 
propriate. The present work is based on a data 
set covering a broad spectrum of objects, includ- 
ing Galactic and extragalactic GCs, UCDs, young 
massive star clusters, nuclei of dwarf ellipticals 
and pressure-supported galaxies, presented in Sec- 
tion 2. We have used the multivariate statistical 
method of K-Means cluster analysis (presented in 
Section 3), to classify these diverse objects with 
respect to a set of physical parameters. Six ho- 
mogeneous groups have been identified by this ob- 
jective method, they are described in Section 4, 
and their properties have been studied by several 
two-point correlations and regressions (Section 5). 
Finally conclusions have been drawn in Section 6. 
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2. The data 

The present sa mple is compo s ed of 370 objects 
from the paper of iForbes et al. (l2008l) . hereafter 
F08, to which we added 19 GCs in M31. We did 
not use all objects of F08 because we were not 
able to document all the values of the additional 
parameters (age, metallicity, colors) used in the 
present study. We took the distance, central ve- 
locity dispersion CTq, effective radius and ap- 
parent K magnitude rriK from Table 1 of F08. We 
did not use the Rf^ values of the galaxies given in 
Table 1 of FOB as they did not agree with Fig. 3 
of FOB. Instead we recalculated R^ using i?20 fol- 
lowing the method outlined in FOB. To that end 
we needed the axis ratio of the galaxies. We ob- 
tained i?25 from the hyperleda databaseQ. logR25 
is the logarithm of the axis ratio at the isophote 25 
mag/arcsec^ in the B band. It was available for 
all galaxies except two (NGC1273 and PER195), 
which were removed from the sample. The Rh 
values that we obtained agree qualitatively with 
those plotted in Fig.3 of FOB (see Fig. [1]). To that 
sample we added 19 GCs in M31. For these addi- 
tion al GCs the stru ctural parameters were taken 
from Peacock et al.l ( 2009) , the velocity dispersion 
from 
from 



Strader et al.l (120 09') and the K magnitude 
Galleti et alT 120041 ) . Hereafter we use the 



term IMO to designate dSphs and what FOB call 
intermediate-mass objects, which include UCDs, 
young massive stellar clusters, nuclei of dEs and 
M32. 

We then derived the virial mass (Myir) us- 
ing the method outlined in FOB, the absolute K 
magnitude Mk and the virial mass-to-light ratio 
in the K band {Myir/Lx)- We derived the ef- 
fective luminosity density 1^ in the K band (in 
Lk,q/pc^) and the effective surface brightness 
^■h,K (in mag / arcsec^) using the relations 



Logle = QA{Mk,q- MK)~log{2TT)~2logRh (1) 



lih^K = MK,Q + 21.b72-2MogIe, (2) 

where Mk,q = 3.2B. 

Next we derived or collected from the literature 
the metallicity, broad-band colors, stellar mass-to- 
light ratio and age. The age is the most poorly 



known parameter of all, except perhaps for the 
Galactic GCs, whose relative ages are well known 
from studies of deep color-niagnitude diagrams 



e.g. iDe Angeli et al.l 120051 : iMarin-Franch et al.l 
20091) . The term "age" is defined precisely only 
for stars and globular clusters, which originated in 
a single star forming burst. For galaxies different 
methods give different age estimates. Integrated 
characteristics, like colors, or narrow-band indices 
are simultaneously influenced by metallicity and 
age, resulting in a degeneracy. In this paper the 
term "age" designates the age of the main star 
formation period, and the metallicity of a galaxy 
is its mean metallicity. Table 1 lists all the pa- 
rameters considered in the present work. 

2.1. Galactic GCs 



Metallicities in the IZinn fc West! (Il984f) scale 
were extracted fro m the McMaster cata log ( Harris! 
20031 ). Ages in the lZinn fc West! (|l984[ ) scale were 
computed from the corre sponding relati ve ages in 
that scale extracted fromlDe Angeli et a l. (200 51), 
Marin-Franch et al.l (|2009l ). lForbes fc Bridgesl(l201^ 
M/Ly were computed using the lBruzual fc Chariot 
(|2003l ) models with the Padova (1994) tracks and 
the Chabrier IMF and GC colors, corrected for 
Galactic extinction. 

2.2. GCs in M31 



Me tallicities wer e (i) e xtracted from the cata- 
log of Galleti et al. ( 20091 ) and (ii) calc ulated us- 
i ng a f ull s et of Lick indices published bv lPuzia"al 



(l2005l) and lBeaslev et al.l(l2004L[2005l) and the pro- 
gr am of interpola t ion an d chi- square minimization 
of Sharina et al. (l2006l) and ISharina fc DavoustI 
(|2009l ) . Ages were calculated using the f ull set 



of L ick indices p u blished by 



and iBeaslev et al.l (|2004 120051) and the pro- 



Puzia al.l (|2005h 



gr am of interpola t ion an d chi-squa re minimization 
of ISharina et al.l (|2006l) and iSharina fc Davousti 
([20091). For clusters without Lick indices we 
used data from SED and fitting of IWang et al.l 
(l2010[ ). Fo r the other clusters only broad-band 



■"-http: / /Icda. univ-lyonl.fr 



colors fromlGalleti et al.l (|2004|) and M g2. Mgb. 
Fe5270, and Fe5335 from lGalleti et al.l (|2009[ ) are 
available. We obtained approximate ages us- 
ing simple stellar population (SSP) models and 
the color/index data. The latter ages are the 
least accurate. M/Ly were estimated using the 
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Bruzual fc Chariot ( 2003 ) model dependence be- 
tween age and M/Ly at a given age and [Fe/H]. 
All the derived data ar e in agreement with the 
parameters published bv lCaldweh et alJ (l201lh . 



2.3. GCs in NGC5128 



Monachesi et al.ll201l[ ) {[Fe/H] = -0.2 dex). The 
differ ence between the early (jPavidge fc Jonesi 
1992h {[Fe/H] = -l.l± 0.2) and later estimates 
is huge. However, a large spread of metallicities 
and ages for stars has recently been found i n M32 



(|Coelho et al.ll2010l: iMonachesi et ai]|201l[ ). The 



T,, . 1,. ■,. i 1 f 1 1 . Ar\ u • Z TOHFl were taken equal to 13 Gyr for all dSphs, 

Metallicities were taken irom lable 4 oi lDabrmghausen et°ar.l ^ •' ^ ' 

because all of them show ancient periods of star 



(|2008[ ). For GCs without metallicity from the lit- 
eratu re we calculated metallic ities using a relation 
from lSalaris fc Cassisil ()2007l) : [Fe/H] = (3.87 ± 
0.07)(y-/)- (5.14±0.08). Ages and M/Ly were 
estimated using broad-band colors from hyper- 
leda, extinction data and the iBruzual fc Chariot 
(|2003l) models. Internal extinction in NGC5128 
was included. 

2.4. IMOs 

Age and metallicities of the young star clusters 
in the remnant of the "wet" merger NGC34 (W3, 
W30, G114) were taken from ISchweizer fc Seitzer 
(|2007[ ). M/Lv were calculated using SSP mod- 
els and photometric results presented in that pa- 
per. We used broad-band photometry, metallic- 
ities, Lick indices and M/Ly for VUCD 3, 4, 5 
from lEvstigneeva et al.l (|2007l ) to calculate their 
approximate age. Evolutionary parameters and 
photometry results of the dwarf-globular transi- 
tion objects in the Virgo c luster (H8005, S 314 t o 



S999) were adopted from iHasegan et al.l (|2005[ ). 

Photometry and metallicity for the four UCDs in 
the Fornax clus t er (U CD2 to 5) were provided 



bv iMieske et al.l (|2006r ). Ages, metallicities and 
M/Lv of the UCDs in the Fornax cluster (F- 
5 to F 59) were determined by IChilingarian et aL 
(|201l[ ). The evolutionary parameters of the UCD 
M59cO were taken from IChilingarian fc MamonI 
12008,). Th e metallicity for the d E VCC1254N was 



taken from lDurrell et al.l (|1996( ). 



Colors and metallicities of dSph s in t he Lo- 
cal Group were taken from iMateol ( 19981 ). The 
peculiar dwarf galaxy M32 was included in the 
list of IMOs (not Es) by F08 with morpholog- 
ical type c E (compact e lliptical). We used the 
data from iMateol (|1998[ ) for all its parameters, 
except metallicity. Since M32 is seen through 
the disk of M31, its very crowded surroundings 
make it a complex case for photometry and spec- 
troscopy. We used the result s of deep CMP studies 
for the metallicity estimate ( Grillmair et al. 19961: 



formation according to CMD studies. 

2.5. Giant and dwarf elliptical galaxies 



105 galaxies from the literature (1 Jerien et al. 


2004; 


Proctor et al.l 


1994 Sanchez-BlazQuez et al. 


2006; 


Li et al. 2007 


: ISerra et al.l 20081 Annibah et al. 


2007: Chilingarianl 20091). These all are 


spec- 



troscopic determinations, except in the first pa- 
per. For many of the sample ellipticals we 
have only broad-band colors corrected for extinc- 
tion from HyperLeda. W e used SSP models of 
Bruzual fc Chariot! (|2003[ ) and a Chabrier IMF to 



derive M/Ly, ages and [Fe/H] using broad- band 
integrated colors and magnitudes. 

Integrated colors not only depend on both 
mean metallicity and age of a galaxy, but may 
also be affected by internal extinction, possi- 
ble ionized gas emission near the galactic cen- 
ter, etc. There is a large number of unknown 
parameters. To model the influence of age and 
metallicity on integrated colors, we use the fun- 
damental luminosity-metallicity re lation common 
for dwarf and gian t ellip ticals (jPrugniel et al 
19931: iThomas et aP l2003l). We used the re- 



sults of simulations of a galaxy with exponen- 
tially declining star burst to derive the approx- 
imate dependence of the broad-band colors on 



M/Ly, age, and [Fe/ff] (jBeh fc de Jongl 12001 
Matkovic fc Guzmanl 20051 ). We selected colors 
more sensitive to age or to metallicity. I — K 
shows a minimal dependence on age and M/Ly. 
There is a strong correlation between B — R and 
stellar M/L ratio i ndependent of metal licity or 
star formation rate ( Bell fc de Jongl 2001 ) . U — B 
is very sensitive to the age of a stellar system and 
its M/L ratio. The slope of the color-magnitude 
relation, and the color - velocity dispersion (cr) 
relation mainly depend on metallicity. Since mass 
is correlated with a for Es, the color - veloc- 
ity dispersion (ct) relation is equivalent to the 
mass metallicity relation. The age difference be- 
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tween galaxies contributes mainly to the scatter 
of the mass-metallicity relation. Fig. 2 compares 
our metallicity estimates for 105 ellipticals of our 
sample with values from the literature. 

Table 1 summarizes the data described in Sec. 2. 
The successive columns give : name, logRh (Rh 
in pc), Mk, fJ.h,K, PUir/LK, M/Lv, the broad- 
band colors (U - B), (B -V), (V - I) and (B - K), 
the metalHcity {[Fe/H]) and age determined by 
us, the metallicity and age from the literature, 
the reference to the latter two data, and finally 
the group. Our main contribution to the data is 
to have derived ages and metallicities for 26 GCs 
in M31, ages and metallicities for Es, and stellar 
mass-to- light ratios for most objects of the sample. 

3. The K-Means clustering technique 

Cluster analysis (CA) is the art of finding 
groups in data. Over the last forty years different 
algorithms and softwares have been developed for 
CA. The choice of a clustering algorithm depends 
both on the type of data available and on the 
particular purpose. 

In the present study we have used the K-Means 
partitioning algorithm (jMacOueenlll967l ) for clus- 
tering. This algorithm constructs K clusters i.e. 
it classifies the data into K groups which together 
satisfy the requirement of a partition such that 
each group must contain at least one object and 
each object must belong to exactly one group. So 
there are at most as many groups as there are ob- 
jects {K < n). Two different clusters cannot have 
an object in common and the K groups together 
add up to the full data set. Partitioning methods 
are applied if one wants to classify the objects into 
K clusters where K is fixed (which should be se- 
lected optimally). The aim is usually to uncover 
a structure that is already present in the data. 
K-Means is probably the most widely applied par- 
titioning clustering technique. 

To perform K-Means clustering we used the 
MINITAB package. The K-means clustering tech- 
nique depends on the choice of initial cluster cen- 
ters. But this effect can be minimized if one 
chooses t he cluster cent ers through group average 
method ( Milligan 1980[ ). As a result, the forma- 
tion of the final groups will not depend heavily on 
the initial choice and hence will remain almost the 
same according to physical properties irrespective 



of initial centers. 

With this algorithm we first determine the 
structures of sub populations (clusters) for varying 
numbers of clusters taking K = 2, 3, 4, etc. Then 
using the Gap Statistics (see below) we determine 
the optimum number of groups. 

3.1. The Gap Statistics 

In order to find the optimum number of 
grou ps we follow the algo rithm of Gap Statis- 
tics ( Tibshirani et al.l 200ll) . Suppose that a data 
set Da, i = 1, 2, n, 1 = 1, 2, p, consists 
of p features measured on n independent observa- 
tions. Let dij denote the distance between obser- 
vations i and j. The squared Euclidean distance 
^liUii — VjiY is used as a most common choice 
for dij . Suppose that the data have been grouped 
into k groups Gl, G2, Gk, with Gr denoting 
the indices of observations in group r, and rir is 
the number of observations in group r. Let 



Dr — ^ ^ dij 



(3) 



be the sum of the pairwise distances for all points 
in cluster r, and let 



(4) 



In the case that d is the squared Euclidean dis- 
tance, Wk will be the pooled within-cluster sum 
of squares. The graph of log{W^) is standardized 
by comparing it with its expectation under an ap- 
propriate null reference distribution of the data. 
The estimate of the optimal number of clusters is 
then the value of k for which log{Wk) falls the far- 
thest below this reference curve. Hence the gap is 
defined by 



Gapnik) = E*Jog{Wk) - logiWk) 



(5) 



where denotes the expectation from the ref- 
erence distribution. The estimate k will be the 
value maximizing Gapn (k) on the basis of the cor- 
responding sampling distribution. As a motiva- 
tion for the Gap Statistics, one may consider clus- 
tering n uniform data points in p dimensions, with 
k centers. Then assuming that the centers align 
themselves in an equally spaced fashion, the ex- 
pectation of log{Wk) is approximately 
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log{pn/12) — {2/p)log{k) + constant (6) 

In other words, the Gap Statistics is defined 
as the difference between the log of the Residual 
Orthogonal Sum of Squared Distances (denoted 
log{Wk)) and its expected value derived using 
bootstrapping under the null hypothesis that there 
is only one cluster. In this implementation, the 
reference distribution used for the bootstrapping 
is a random uniform hypercube, transformed by 
the principal components of the underlying data 
set. If the data actually have K well-separated 
clusters, then it is expected that log{Wk) will de- 
crease faster than its expected rate (2/p)log(k) for 
k < K. When k )^ K, then a cluster center 
is essentially added in the middle of an approx- 
imately uniform cloud and simple algebra shows 
that log{Wk) should decrease more slowly than its 
expected rate. Hence the Gap Statistics should be 
largest when k = K. 

3.2. The algorithm to find the Gap Statis- 
tics 

Two common choices for the reference distri- 
bution are : (a) each reference feature is gener- 
ated uniformly over the range of the observed val- 
ues for that feature; (b) the reference features are 
generated from a uniform distribution over a box 
aligned with the principal components of the data. 

In other words, if X is an nxp data matrix, 
it is assumed that the columns have mean and 
then the singular value decomposition X = UD^-^ 
is performed. It is transformed through Y = XV 
and then uniform features, say T, are drawn over 
the ranges of the columns of Y, as in method (a) 
above. Finally it is back-transformed via Z = TV^ 
to give reference data, say Z. Method (a) has the 
advantage of simplicity. Method (b) takes into 
account the shape of the data distribution and 
makes rotationally invariant, as long as the clus- 
tering method itself is invariant. 

In each case, E^log{Wk) is estimated by an av- 
erage of B copies log{Wj^), each of which is com- 
puted from a Monte Garlo sample Y*,Y2, ■■■,¥* 
drawn from the chosen reference distribution. Fi- 
nally, one needs to access the sampling distribu- 
tion of the Gap Statistics. Let Sd{k) denote the 
standard deviation of the B Monte Carlo replicates 



log{W^). Accounting additionally for the simula- 
tion error in E^log{Wk) results in the quantity 
Sfc = + l/B)sd{k). Using this the estimated 
cluster size k is chosen to be the smallest k such 
that Gap(k) > Gap(k + 1) — Sk+i- where Sk+i is a 
function of standard deviation of the bootstrapped 
estimates. 

The computation of the Gap Statistics proceeds 
as follows: 

• Step 1: The observed data is clustered by 
varying the total number of clusters from 
k= 1, 2, K, giving within-dispersion mea- 
sures Wk, k = 1, 2, K. 

• Step 2: B reference data sets are gener- 
ated using the uniform prescription (a) or 
(b) above and each one is clustered giving 
within-dispersion measures W^^^, b = 1, 2, 

B, k = 1, 2, K. Then the estimated Gap 
Statistics is calculated as follows: Gap(k) = 
{l/B)Y.,log{W;,)-log{Wk). 

• Step 3. Let /"= Y.h '0.9(^^4), then the 
standard deviation is computed as 

sdkH{l/B)Y.,log{W*^,)-f]^ 

and Sk is defined as Sk=sdk\J {I + 

Finally that number of clusters are chosen such 
that K;=smallest k and 

Gapik) > Gapik + 1) - Sk+i (7) 

In other words the optimum number of clusters 
is that k for which the difference 



Uk = Gap{k) - {Gapik + 1)- Sk+i) > (8) 

4. Results 

The parameter set chosen for CA consists 
of Mk, log{ao), logRu, Myir/Lx, [Fe/H\,M/Lv. 
Parameters like Myir, iJ,h,K and age are not used. 

Mk 

is very highly correlated with Myir and iih,K, 
so inclusion of these parameters docs not influ- 
ence the clustering. Age is excluded because of the 
large uncertainties associated to it. The remain- 
ing parameters are not used because of the large 
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number of missing values. But, once the substruc- 
tures are identified, all the parameters are used to 
identify the distinctive properties of the groups. 

We have calculated the Gap Statistics for the 
set of six above parameters, and the output sug- 
gests that the optimum number of clusters is ei- 
ther four or six because the criterion used in Gap 
Statistics to find the optimal number of clusters, 
i.e. Uk — Gap{k) — {Gap{k -I- 1) — Sk+i) > is sat- 
isfied for k = 4 and k = 6. Table 2 and Fig. 3 show 
that the value of exceeds for k = 4 and k = 6 
and there is a sharp decline of the graph after the 
value k = 6. Hence, considering all the criteria 
discussed above, the optimal number of clusters 
for the present sample is k = K = 6. The six clus- 
ters (hereafter named groups to avoid confusion 
with star clusters) are designated FKl to FK6, 
and their average properties are given in Table 2. 

The elliptical galaxies were divided into two 
groups by the CA: high-mass ellipticals (gEs) in 
FKl and low-mass ones (dEs) in FK5. Note that 
the labels "gE" and "dE" do not refer strictly to 
the morphological types commonly used in astron- 
omy. We use these designations conditionally, to 
stress the statistical difference in mass between the 
objects of FKl and FK5. 

FK3 has the high-metallicity GCs and the 
bright and high-metallicity IMOs. The bright- 
est IMOs are UCD2 [Mk = -16.32), VUCD3 
[Mk = -16.21) and UCD3 (Mk = -1 6.215) 
Iffilker et al.ll2007t lEvstigneeva et al.ll2007t) . 

FK6 is composed of IMOs, of the most mas- 
sive GCs in the Galaxy, in M31, and in NGC5128, 
(these GCs are all of low metallicity), and of dSphs 
of the Local group : Leo I (Dist.= 0.25 Mpc) and 
Sculptor (Dist.= 0.08 Mpc). 

Two groups, FK2 and FK4, have a negligibly 
small number of members compared to the other 
groups : they contain three members each. These 
are Local Group dSphs, listed according to their 
distances from the Sun in Mpc: UMi (FK2, 0.066), 
Draco (FK2, 0.086), Sextans (FK2, 0.086), Ca- 
rina (FK4, 0.1), Fornax (FK4, 0.14), Leo II (FK4, 
0.21). We unfortunately do not have the full set of 
parameters for the other dSphs in the Local Group 
and nearby groups to include them in the analy- 
sis. A probable reason why these six objects were 
classified in such a way is their M/L ratio, which 
is higher than for the other galaxies. These two 



groups are considered only briefly, as their study is 
the subject of a separate and elaborate study. So 
there are essentially four groups found as a result 
of our cluster analysis. 

We show in Fig. 4 how the various types of ob- 
jects are distributed among the different groups in 
Mk - logRh space. The six groups are indicated 
by different symbols, colored according to the mor- 
phological type of the objects: GCs in black, IMOs 
in green, and ellipticals in red. 

To justify our choice of sample, and to show 
that it is representative of stellar systems in 
the local Univers e , we compare it to that of 
(|Misgeld k Hilkeil boill ). (hereafter MH2011) 



who, like us and F08, studied a sample of stel- 
lar systems covering a large range in masses, sizes 
and luminosities. 

The sample of MH2011 is larger than ours, but 
the associated data do not include velocity dis- 
persions, metallicities or ages, so we could not 
perform a similar analysis with their data. Nev- 
ertheless, our sample covers basically the same 
space in absolute magnitudes and effective radii, 
as shown in Fig. 4, which can be compared to Fig.l 
of MH2011. The main difference is that the sam- 
ple of MII2011 has many more dEs (in the Hy- 
dra I and Centaurus clusters of galaxies), and ex- 
tragalactic GCs (mostly GC candidates in Virgo) 
and they included much fainter dwarf galaxies of 
the Local Group, for which velocity dispersions 
would be very difficult to measure. In short, our 
sample does not appear to be biased against any 
particular type of object. 

We also computed the probability density dis- 
tribution (PDF) of Mk in our sample and com- 
pared it to the same distribution for the MH2011 
sample (see Fig. 5). The method of non-parametric 
density estimates is describe d in a previous paper 
(|Chattopadhvav et al.ll2009l ). The bin width for 
computing the density estimates is the same as 
for the histograms shown in Fig. 5. Since MH2011 
gives Mv rather than Mk, we simply shifted their 
V magnitudes by 2.90, which is the average value 
of (V - K) in our sample. There are three main 
populations in both samples, the faintest one be- 
ing much more important in MH2011. Antici- 
pating on our results, we expect the distribution 
of metallicities for the F08 and MH2011 sam- 
ples to be similar due to the fact that Es fol- 
low the fundamental luminosity-metallicity rela- 
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tion (jPrugniel et al.lll993HThomas et al.ll2003l ). 

The successive peaks are at Mk = -24.8, - 
19.7, -13.3, -11 in our sample, and at Mk = - 
25, -19.4, -14.8 in the sample of MH2011. The 
first peak in both samples corresponds to bulges 
and the brightest elliptical galaxies. The next 
peak appears at the location where the linear size 
- luminosity relation, common for ellipticals and 
UCDs (MH2011), splits into two : one relation 
for dwarf galaxies and one for compact ellipti- 
cals and GCs. This occurs at about Rh = 1.3 
kpc and Mass — 10^'^ Mq. So, galaxies in this 
group have roughly constant effective radii. The 
faintest objects in this group have luminosities 
similar to M32, Mk ~ 18.5, but their stellar den- 
sities are two orders of magnitude lower (see Fig. 5 
in MH2011). The highest stehar density for this 
group may be a characteristic scale, dividing stel- 
lar systems into two systems. The internal acceler- 
ation for one group is within the limits postulated 
in MONDian dynamics, while for the other groups 
it is outside those limits. See also the caption of 
Fig. 7 of MH2011. The faintest broad PDF peaks 
(-13.3 and -11 in our sample and -14.8 in MH2011) 
are different for both samples. However, this is 
just a selection effect : as mentioned above, our 
sample contains fewer dEs. 

So, again, our sample does not differ signifi- 
cantly from another large sample of stellar sys- 
tems. Our sample does not reflect the local lu- 
minosity function for individual types of objects, 
and neither does the sample of MH2011. We sug- 
gest that the relative intensity of the PDF peaks in 
both samples reflects the way in which the samples 
were selected. 

We also examined whether our choice of objects 
in the F08 sample (370 out of 499) could bias the 
results in some way. We have computed the mean 
± standard error values of Mk and log{ao) in the 
sub samples 1,2, and 3 considered by F08 as well 
as for our corresponding sub samples. The num- 
ber of objects is of course different in the present 
sample and in the F08 sample. But from Table 4 
it is quite clear that this feature does not intro- 
duce any significant bias as the mean values are 
very similar. 

We now present the distinctive properties of 
the groups, and look for possible physical rea- 
sons for the differences and similarities between 
the groups. 



5. Properties of the groups 

5.1. Mass-to-luminosity ratios and bind- 
ing energies 

We will discuss the virial M/L ratio, and it is 
important for what follows to keep in mind that 
the stellar M/Ly derived using photometric data 
and SSP models is not necessarily identical to the 
true baryonic M/L. This is due to the difficulty 
to correctly take into account the star formation 
history (SFH) and initial mass functi on of stel- 
lar populations (e.g. lTrager et al.ll2008l . MH2011). 
Furthermore, a disagreement between virial and 
baryonic M/L may be due to the presence of dark 
matter, if the stellar population model including 
SFH and initial mass function is correct. 

The difference between virial and stellar M/L 
for our sample can be seen from Table 1. It is seen 
that both the virial {Myir / L k) and the baryonic 
(M/Lv) mass-to-light ratios differ at a high level 
of significance among the four main groups. Here- 
after we will concentrate on Myir/LK and simply 
call it M/L. It is well know n that UCDs tend to 
have higher M/L than GCs ( Dabringhausen et al.l 
20081 . F08), and that dwarf spheroidal galaxies 
have very high M/L from direct radial veloc- 
ity measurements o f their brightest stars (e.g. 
Simon fc Gehall2n07h . Additionally, UCDs, like 
galaxies, have relaxation tim es greater than the 
Hubble time ([Kroupa Il998l) . This is usually 
demonstrated by pl otting the data in th e ki — 
space introduced bylBenderelaD (|l992h . and this 
is well discussed in the aforementioned papers. 

For the present data, these parameters are : 

fci = {logal + logRh)/V2 
k2 — {logaQ + 2logIe — logRh) / V& 



and 



^3 = [logal - logle - logRh)/Vi 



where Je is given by Eq. 1. These coordinates are 
simply related to physical quantities : ki is propor- 
tional to the logarithm of mass, k2 is proportional 
to the effective surface brightness times M/L, and 
k^ is proportional to the logarithm of M/L. 

The differences in mass (represented by fci ) and 
M/L (represented by k^) between the groups are 



8 



shown in Fig. 6. The groups occupy different loca- 
tions in this projection of the FP, except FK3 and 
FK6. For these two groups there is no continuity 
break in the fci, fca parameter distributions as for 
other groups. Both FK3 and FK6 contain objects 
with high M/L. FK3 includes IMOs, and FK6 con- 
tains dSphs (Sculptor and Leo I) and IMOs. We 
also note that the four main groups show wide and 
different ranges in both mass and M/L. In each 
group, more massive objects show higher M/L, 
but the slope of the correlation is different for each 
group. 

To quantify this, we performed robust multilin- 
ear regressions of the form = a + bki on the 
four main groups. The resulting fits are listed 
in Table 5. The regression lines for the groups 
FK3, FK5, and FK6 correspond to the relation 
M/Loc M°-2 within the errors. M/L is propor- 
tional to M° '^^ for the group FKl. The position 
of the different objects within the groups on the 
FP reflects not only differences in M/L, but also in 
surfa ce density, luminosity, and kinematical struc- 
ture ( Djorgovski fc Davis 19871 ). According to the 
slopes of the relations, the objects in FKl are 
much more influenced by the above three factors 
than the objects in FK3, FK5 and FK6. 

We now move on to discuss the edge-on projec- 
tion o f the Fundamental Plane ( Diorgovski fc David 
19871: iFaber fc JacksonI Il976l : iKormendvl 119771 
Diorgovskil Il995( ) shown in Fig. 7. This fig- 
ure is a representation of the Virial Theorem: 
usually applied to g alaxies 



r. DC a?,l7HM/L) 



(|Faber et al.Ml989l : iDiorgovski et al.lll989f l. This 
figure al so serves to compar e the binding energies 
of GCs (lMcLaughlinll2000l ). The most compact 
and luminous GCs have larger binding energies. 

The groups FK3, FK5, and FK6 (i.e. GCs, 
IMOs and dEs) follow roughly the same relation 
in the edge-on projection of the FP (Fig. 7). We 
obtained a bivariate least squares solution fitted 
through fih,K- 

logRh-2logao = 0.4(1. 07(±0.03)/i,i,/f-f 19. 1(±0.1)) 

which corresponds to oc ctqI'^-^. The gEs 
of FKl are concentrated in a parallel sequence, 
shifted towards lower surface brightnesses (nh.K)- 
The bivariate correlation for FKl gives: 

logRh~2logcjQ = 0.4(1. 04(±0.06)/i/i,K+20.0(±0. 13)) 



These two solutions are close to the one that 
satisfies the Virial Theorem. The different slope 
(1.07 in the first case, 1.04 in the second) is re- 
ferred to as the tilt in the FP, whose cause i s still 
under debate (see iFraix-Burnet et al l2010l and 
references therein). The tilt of the virial mass - to- 
tal stellar mass relation common for gEs, cEs and 
UCDs/GCs has been discussed in F08. The differ- 
ence in__tlie_zero points includes three components 
(e.g. lKormendvl[l989l and references therein) . The 
first one refiects the density, luminosity and kine- 
matic structure of objects. The second factor indi- 
cates whether the system is gravitationally bound 
or virialized. If the deviation from the FP is due 
to mass-to-light ratio, this implies the scaling rela- 
tion M/L (X Af°-^. The systematic shift between 
gEs and the groups of GCs, IMOs, and dEs is 
mainly due t o the approximately ten tim es larger 
M/L for gEs (jDabringhausen et aLlliooi l. 

The objects of FK3 and FK6 are well mixed to- 
gether in Fig. 7, with a tendency for FK3, which 
contains IMOs, to have higher binding energy. 
The objects with the strongest deviation from the 
relation are IMOs: e.g. BOOl, M59cO, UCD3; 
the globular cluster NGC2419, and some dEs, like 
IC3779, with ^ih.K > 20 mag arcsec"^. M32 has a 
very high binding energy, similar to that of IMOs. 
Some gEs also fall in the same regio n of the dia 



gram as M 32, but no other dE does. iBekki et al 



([20011) and I Graham! (120021 ) argued that M32 is 
the stripped core of a larger galaxy. NCG2419 
shows a lower binding ene rgy than other GCs. 

considered it as the 



iDabringhau sen et al.l ( |200S 
most likely candidate to host dark matter. 

The ^J.h,K versus logRh diagram (Fig. 8) illus- 
trates the difference in stellar densities between 
the GCs of FK3 and FK6. It shows that the GCs 
in FK3 have higher hh.k than those in FK6 at 
a given Rh- In other words, FK6 has statistically 
shallower surface brightness profiles than FK3. On 
the other hand, the IMOs and GCs in FK3 are 
more massive/lumino us and compact in general 



than those in FK6. [Jordan et al.l (|2005l ) found 



a significant correlation between half-light radius 
and color for early-type galaxies in the Virgo clus- 
ter in the sense that the red GCs are smaller than 
the blue ones. 

Having studied how mass is related to luminos- 
ity in our different groups, we now examine how 
mass is related to metallicity. 
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5.2. Mass-metallicity relation 

It is now well established that more massive 
galaxies are also more metal rich; this is a conse- 
quence of the hierarchical formation of galaxies in 
the Universe. But does such a relation hold for all 
types of stellar systems? 

5.2.1. A boundary line 

The mass-metallicity relation (hereafter MMR) 
for our sample objects is shown in Fig. 9a, where 
ki, which is equivalent to mass, is plotted versus 
[Fe/H]. 

This figure shows that, except for a few objects, 
all types of stellar systems lie above a boundary 
line. It was plotted to stress the tendency, but its 
slope is surprisingly close to a MMR of the form 
Z cx Mass'^'^. The correlation is very weak for the 
objects in FK3 {r{MK,[Fe/H]) = -0.382) and 
FK6 ir{MK, [Fe/H]) = -0.177), if we consider 
them globally. Only the brightest low-metallicity 
GCs (FK6) and IMOs at a given metallicity are 
close to the MMR. The picture is almost the same 
if we plot absolute K magnitude versus [Fe/H] 
(Fig. 9b). However, here the slope of the boundary 
line is slightly different from that of the MMR: 
[Z/H] ~ -3.5 - O.UMk- 

What could be the origin of the boundary 
line? It is unlikely to be caused by an obser- 
vational selection effect. We would presumably 
not see it if we included in the sample only high- 
metallicity GCs and galaxies of other morpholog- 
ical types. Many of the GCs are the brightest 
GCs of our Galaxy, and their metallicities are very 
accurate. Extragalactic GCs and IMOs are also 
bright. Their metallicities were obtained mainly 
via spectroscopy, and are not very much influenced 
by observational errors and the age-metallicity de- 
generacy. The only really uncertain metallici- 
ties are those of ellipticals, because of the age- 
metallicity degeneracy and uncertainties due to 
possible internal extinction, light-element abun- 
dance variations, and large age and metallicity 
spreads within individual galaxies. But, in spite 
of these uncertainties, the ellipticals do follow the 
relation. 

The slope of the boundary line is similar to 
that of the luminosity-metallicity relation found 
in the literature. A luminosity-metallicity rela- 
tion, [Z/H] ^ -3.6 - 0.19Mb, was found for 



dwarf a nd giant ellipt i cals in nearby galaxy clus- 
ters by I Thomas et al. (l2003h . It is equivalent to 
the equation Z oc found for dw a rf gal axies in 
the Local Group by iDekel fc Silk (Il986), since 
\Z/H ] = [Fe/H] + 0.94[a/f e] (iThomas et all 
l2003h . log{LB/LQ) = .4(5.48 -Mr) , and I oqZ • 



0.977[Fe/H] - 1.699 (jBertelh et al.lll994l) . We 



used here the solar value [a/Fe] = 0. However, 
the deviations from this relation for massive el- 
lipticals may be large due to strong variations in 
[a/Fe] for Es: ^ .2-^0.5 dex (|Thomas et al.ll2003 : 
Puzia et al.l l2006l and references therein) . The 



median of the metallicity distribution for elliptical 
galaxies and galactic bu lges from the Sloan El igital 
Sky Survey obtained bv lGallazzi et al.l ()2005h as a 
function of stellar m ass is also close to the relati on 
Z cx M" "* (see also iDabringhausen et al.|[2008l) . 

The origin of the luminosity-metallicity and 
mass-metallicity relations for different morpho- 
logic al types of galaxies is still an open issue 



Grebel et al, 2003: ,Finlator fc Dave I l2008t 



Kunth fc Ostlinll2000l) . Does star formation de- 



fine the shape of the MMR? Does the boundary 
line mean a lower fraction of matter capable of be- 
ing transformed into stars under special physical 
conditions? It might result from the interplay be- 
tween internal and environmental factors: mergers 
and interactions, inflows and outflows of gas, star 
formation histories of individual galaxies in hier- 
archical galaxy formation. 

The luminosity-metallicity relatio n for bright 



est GCs has been extensively studied [Harri s et al 



20061: iMieske et all 120061: iPeng et al., ,2009il . Us 



ing linear color-metallicity relations for blue GCs, 
these studies derive scaling relations between GC 
luminosity L and metallicity Z con sistent with 
Z cx fe.g. IStrader fc SmithI lioosi) . The slope 
depends on the SSP models and on the light- 
element abundances. According to Carnev ( 1996h 
the mean [a/Fe] for Gala ctic GCs is 0.3 dex 
The s ame value was used by IDabringhausen et al 



(|2008[ ) to calculate [Z/H] for IMOs. 

The metallicity of the faintest GCs close to the 
MMR is intriguing. It corresponds approximately 
to extreme abundances of Population II stars, i.e. 
stars formed immediately after the initial pollu- 
tion of interstellar med ium by ma ssive Population 
HI stars: Z - O.OIZq (|Silklll985l) . 

The Color-magnitude diagram (CMD) and 
chemical composition of some GCs located near 
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the border line (i.e. flCen, NGC2419, NGC 6341) 
are unusual. For example, NGC2419 is consid- 
ered a remnant of a dwarf g alaxy due to its pe - 



culiar chemical composition (jCohen et al.ll201Cl() . 
The CMDs of most of these GCs show the exis- 
tence of multiple stellar populations, a fact that is 
still n ot fully understood (see e. g. iBedin et al. 
20081 : iMarin-Franch et all l2009l : iRenzinil l2o'o8l 
and references therein). Since these GCs and 
IMOs are close to the MMR, they were probably 
the brightest parts (nuclei) of tidally destructed 
host galaxies ( Zinnecker et al.l Il98 1 iBekki et al.l 
20031 ). They may also be geniune compact 
dwarf galaxies originating from small-scale peaks 
in the primordial dark matter power spectrum 
([Drinkwater et al. jCs may have formed i n 

dark matter minihaloes ( Mashchenko et al. 2005[) . 
However, it has not been established whether they 



actually contain dark rnatter haloes (|Jordi et al. 
20091: iBaumgardt et aklboOOl ). 



The Local Group dSphs and some GCs defi- 
nitely fall below the border line in the fcl versus 
[Fe/H] diagram. The reason has been studied ex- 
tensively for dSphs. Dwarf galaxies with luminosi- 
ties below some limit lose gas effectively because 
of their low gravitational potentials, too shallow 
to prevent stellar outflows following star forma- 
tion e pisodes (e.g. Dekel fc Silklll986 : lGrebel et al 
20031 ). 



5.2.2. Metallicity bimodality 

There is a gap between the groups FK6 and 
FK3 in Figs.9(a, b). It is located near [Fe/H] = 
— 1.0 and is not horizontal. The low- metallicity 
peak is near [Fe/H] = — 1.6 ± 0.4 dex, the high- 
metallicity one is near [Fe/H] = — 0.6± 0.04 dex. 
Similar metallicity peaks and the dividing line 
wer e identi f ied fo r the GC system of our Galaxy 
by ( Harris! Il989l ). GC systems of massive Es 



and many spirals follow a bimodal color distribu- 



tion ( 


Harris et al.lll996; 


Gebhardt & Kissler-PatisJ 


19991 Larsen et al. 2001 


: Peng et al. 20061). Anew 



feature shown in Fig. 9 is that the metallicity dis- 
tributions of IMOs and extragalactic GCs fall in 
the same range as the GCs of our Galaxy. Extra- 
galactic objects are gathered in two homogeneous 
groups together with Galactic GCs. 

Fig. 10 shows the probability density distribu- 
tion of [Fe/H] for all galaxies (FKl -f- FK5 -I- FK4 
-I- FK2), and for GCs and UCDs (FK6 + FK3). 



The distribution is computed in the same way as 
the PDF shown in Fig. 5. Although the groups are 
not plotted with different colors, they are clearly 
distinguishable from the PDF peaks. The local 
maximum in the distribution for group FK3 at 

[Fe/H] 0.6 is close to that of FK5 ("dEs"), 

which is composed of galaxies having a roughly 
constant effective radius and departing from the 
size - mass relation common for gEs and UCDs 
(MH2011). The corresponding PDF peak is also 
present in the probability density distribution of 
luminosities (see discussion at the end of Sec. 4). 
The PDF for group FK6 corresponds to that of 
groups FK2-f 4, i.e. galaxies less massive than 10® 
stellar masses. Interestingly, there is another gap 

at the level of [Fe/H] 0.3, between six IMOs 

(M59cO, W3, W30, G114, VUCD 3, S490) and 
the other objects in FK3. M32 with a central ve- 
locity dispersion of uo ^ 79 km/s also falls in this 
metallicity range. However, since the luminosity 
and metallicity distributions of galaxies are influ- 
enced by sample selection effect, the correlation is 
not sufficient to establish the tidal origin of nuclear 
GCs and UCDs. 

The nature of the bimodality in the metallicity 
distribution is a complex, still unanswered ques- 
tion. Due to the stochastic nature of galaxy for- 
mation and star formation, hierarchical scenar- 
ios do not reproduce the metallicity bimodality 
well. The dependence of the ga lactic SFH on 
stella r mass is not st raightforward ( Thomas et al.l 
20051: iRenzinil l2009t ). Additionally, there is a 
strong morphology-density relation, the environ- 
mental dependence between stellar mass, struc- 
ture, star formation and nuclear a ctivity in galax- 



ies 



(e.g. iKauffmann et al.l 12004 IRenzinil |2006|) 



Recent spectroscopic studies have revealed strong 
age and metallicity gradients of different slopes 
and values between t he nuclear and oute r regions 



of elliptical galaxies (|Koleva et al.ll201lL and ref- 
erences therein). Nuclear activity in galaxies of- 
ten continues longer than in the outer regions, 
because the fuel for star formation falls towards 
the gravitational center. So, nuclei may contain 
multiple stellar populations, and be on average 
younger and more metal-rich than the rest of the 
galaxy. There are anomalous objects in both FK3 
and FK6. In the low-metallicity group GCs like 
nCen, NGC2419 and NGC 6341 show evidence 
of multiple stellar populations. The metal-rich 
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GCs NGC6441 and NGC6338 have pro minent blue 



extensions in the horizontal branch (Rich et al 



19971), which are not typically associated with a 
globular cluster of this metallicity, like 47 Tuc. 

It has been proposed that the outer-halo 
GCs of the Galaxy were accreted from the 
satellite galaxies (e.g. Mackev fc Gilmord 2004t 



from the paper of iForbes et al.l (|2008l ). It consists 
of elliptical galaxies, intermediate mass objects. 
Local Group dwarf spheroidals, nuclei of dwarf el- 
lipticals, young massive objects and globular clus- 
ters. The sample properties were completed by 
data from the literature or derived by us. Our 
aim was to investigate the existence of intercon- 



Chattopadhvav fc Chatto padhvav ll2007l :IChatto mdhvav ei^ dness, if any, among the six groups found by 



Shapiro et al. (201of ) 



2007; .Mondal et al. .2008,) . 
suggest that high-metallicity old GCs were formed 
from super star-forming clumps with radii 1-3 kpc 
and masses 10^ to lO^M©, which are known as a 
key component of star- forming galaxies at z ~ 2. 



5.2.3. Efficiency of metal production 

Is there a similar physical quantity for dynami- 
cally hot stellar systems lying close to the MMR? 
For young stellar systems, for example, the star 
formation rate is known to be an important fac- 



tor influencing the MMR (|Mannucci et al.ll2010l 



and references therein). In Fig. 11 we plot the de- 
pendence of [Fe/H] on the metallicity per unit 
effective luminosity density in the K band. We 
call the last term "metal production efficiency" 
(MPE) by analogy with the star formation effi- 
ciency (SFE), which is the fraction of gas con- 
verted into s tars at a particul ar evolutionary stage 
of galaxies (jKennicuttI 1 19981 ). MPE also reflects 
the stellar density and the size of dynamically 
hot stellar systems. Fig. 11 shows that GCs and 
UCDs in FK6 have MPE in the same range as Es: 
gEs (FKl) (MPE= 2.7 ± 0.25), and dEs (FK5) 
(MPE= 1.7 ± 0.4). Galaxies with stellar masses 
M < 101°Mq, including dSphs and GCs + UCDs 
(FK3-I-FK6), are in two separate sequences, both 
showing a tendency for metallicity to increase lin- 
early with MPE. The two sequences intersect at 
the location of the brightest UCDs and M32-like 
objects. Fig. II also shows that the objects of FK3 
are the most metal-rich per unit effective luminos- 
ity density. So, at least for these GCs and UCDs 
in our sample, it is reasonable to assume that they 
are the densest parts of galaxies accumulating fuel 
for star formation. 

6. Conclusion 

A multivariate statistical technique, K-Means 
clustering, has been carried out on a data set taken 



our multivariate analysis. 

In order to inquire into the physical origin of 
IMOs, we considered different projections of the 
fundamental plane using the results of the statis- 
tical analysis along with observational data on ve- 
locity dispersion, effective radii and effective sur- 
face brightness calculated from the total absolute 
magnitude in the K band. We found that our 
groups are different in terms of virial M/L ratios, 
and dependences between virial M/L ratios and 
mass. 

The value of our study is that we include metal- 
licities along with other data in addition to the list 
of parameters of F08, which definitely helps us to 
provide an objective classification into groups. We 
consider a unified mass-metallicity dependence for 
all the sample objects. It shows that (i) there are 
GCs and UCDs in the low-metallicity group shar- 
ing MMR with galaxies; (ii) there are signatures 
of bimodality/multimodality in the metallicity dis- 
tribution that are common for GCs and IMOs on 
one hand, and for low- and high-mass Es on the 
other hand. We speculate that the rate of SF at 
the epoch when the objects were young is the prob- 
able reason for the above two features. It appears 
that the mean metallicities per effective K-band 
luminosity density (MPE) for GCs and UCDs in 
FK6 lie in the same range as for elliptical galax- 
ies, suggesting similar physical processes and SFE. 
However, MPE is much higher for GCs and UCDs 
in FK3. This confirms that these objects origi- 
nated as the densest parts of the present day Es. 

According to our findings, IMOs may be di- 
vided into two physical groups: (i) Dwarf galaxy 
- globular cluster transition objects formed in the 
same way and from the same material as old galax- 
ies and (ii) nuclei stripped from dwarf and nor- 
mal ellipticals during their dynamical evolution in 
groups and clusters. Note that since UCDs were 
found only in dense environments, the last sugges- 
tion is highly probable. Extensive theoretical and 
observational studies are needed to establish the 
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reasons for the described features and the exact 
nature of UCDs. 
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Fig. 1. — Effective radius Rh computed by us 
(filled circles) and i?20 from F08 (crosses) versus 
absolute K magnitude Mk- This figure is to be 
compared to Fig. 3 of F08. 
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Fig. 2. — Comparison of our derived metallicities 
of Es with values from the literature (see Table [T]) 
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Fig. 4. — Absolute K magnitude versus logarithm 
of effective radius. The elliptical galaxies are in 
red, the IMOs are in green and the GCs are in 
black. Open circles are for FKl, pluses are for 
FK2, crosses are for FK3, asterisks are for FK4, 
triangles are for FK5 and open squares are for 
FK6. 
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Fig. 5. — Probability density functions of absolute 
K magnitude Mk- They are plotted in green for 
our sample and in black for MH2011. The hues 
are non-parametric density approximations. They 
are in red for our sample and in gray for MH2011. 
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Fig. 6. — Projection of the FP for dynamically hot 
stellar systems, fci and are related to mass and 
M/L respectively. The symbols and colors are the 
same as in FiglH 
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Fig. 7. — Edge-on projection of the Fundamental 
Plane for the six groups. The symbols and colors 
are the same as in Fig|4l 
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Fig. 8. — The surface brightness {nh,K) versus half 
hght radius (logRh) profile for GCs in groups FK3 
(green, plus) and FK6 (black, box) respectively. 
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Fig. 9. — ki versus [Fe/H] (a) and Mk versus 
[Fe/H] (b) for the six groups. The symbols and 
colors are the same as in Fig. 4. Black dots for el- 
lipticals indicate metallicities from the literature. 
Black dots for GCs (open squares) indicate Galac- 
tic GCs 
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Fig. 10. — Probability density functions of metal- 
licity [Fe / H] . They are plotted in green for galax- 
ies and in black for GCs and IMOs. The lines 
are non-parametric density approximations. They 
are in blue for galaxies and in brown for GCs and 
IMOs. 
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Fig. 11. — [Fe/H] versus metallicity per unit ef- 
fective luminosity density. The symbols and col- 
ors are the same as in Fig. 9. Bold symbols indi- 
cate that GCs and IMOs tend to follow the mass- 
metallicity relation. 
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Table 1 

Photometric and structural parameters. 
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Table 2 

Gap values for different values of the number k of clusters. 



Number of clusters 


Gap 




k=l 


1.658168 


-0.3570766 


k=2 


1.592197 


-0.1841725 


k=3 


1.804713 


-0.7284219 


k=4 


2.592263 


0.0873069 


k=5 


2.529849 


-0.5264599 


k=6 


2.43953 


0.0961393 


k=7 


2.356531 


0.08300998 


k=8 


2.572305 


-0.4335349 



Table 3 

Average properties with standard errors of the four main groups. 



Groups 



FKl 



FK3 



FK5 



FK6 



Number of members 
logaQ{kms^^) 
logRhipc) 
MK(mag) 

tJ'h,K i'tnag arcsec~^) 

M^ir./LK{M^ir.^Q/LK,e) 

M/Lv(MQ/Lv,e) 

loe(M^ir)iMQ) 

U — B{mag) 

B - V(mag) 

V — I{mag) 

B — K(mag) 

[Fe/H\{dex) 

Age (Gyr) 



210 

2.3385 ± 0.0076 
3.8646± 0.0228 

-24.691 ± 0.058 
18.199 ± 0.059 
3.171 ± 0.116 

I. 9073 ± 0.049 

II. 625 ± 0.034 
0.466 ±0.00541 

0.893 ± 0.003 
1.0001 ± 0.0034 

3.818 ± 0.017 
0.04495 ± 0.0087 

5.904 ± 0.159 



57 

1.3134 ± 0.0329 
0.7363 ± 0.0537 
-13.305 ± 0.214 
13.944 ± 0.169 
1.772 ± 0.239 
3.173 ± 0.290 
6.753 ± 0.106 
0.3205 ± 0.0419 
0.793 ± 0.032 
1.0527 ± 0.027 
3.732 ± 0.270 
-0.5472 ± 0.038 
10.637 ± 0.468 



39 

1.6088 ± 0.0214 
3.013 ± 0.0333 
-19.776 ± 0.143 
18.857 ± 0.176 
2.489 ± 0.421 
1.182 ± 0.170 
9.6175 ± 0.0957 
0.3575 ± 0.0378 
0.870 ± 0.034 
0.9375 ± 0.022 
2.903 ± 0.066 
-0.6069 ± 0.021 
4.817 ± 0.646 



77 

0.9679 ± 0.0346 
0.6397 ± 0.0454 
-11.229 ± 0.198 

15.537 ± 0.177 
1.750 ± 0.147 

2.716 ± 0.230 
5.9656 ± 0.0945 
0.0671 ± 0.0150 
0.689 ± 0.01 

0.889 ± 0.0078 

3.299 ± 0.100 
-1.5030 ± 0.380 

10.217 ± 0.226 



Table 4 



Average properties with stand ard errors of three subsamples taken from IForbes et al. 

( 20081 ) AND from our smaller sample. 



Samples 


M7r(F08) 


Mk (ours) 


(TO (F08) 


(To (Ours) 


1 (GCs) 


-11.327± 1.87 


-11.339±0.156 


12.4±0.841 


12.077±0.686 


2 (IMOs) 


-13.98±0.291 


-14.145±0.246 


28.58±2.45 


29.13 ± 2.44 


3 (Es) 


-23.67 ± 0.084 


-23.943 ±0.124 


185.5 ±3.84 


196.81 ± 5.31 



Table 5 

Robust multivariate regression on the four main groups in kl-ki space. 



Group a b rms 



FKl 
FK3 
FK5 
FK6 



-0.567413 
0.238372 
0.133827 
0.378629 



0.314669 
0.244655 
0.141191 
0.155797 



0.06823 
0.09654 
0.17752 
0.11949 
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